Resonant production of the sterile neutrino dark matter and fine-tunings in the z/MSM 
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The generation of lepton asymmetry below the electroweak scale has a considerable impact on 
production of dark matter sterile neutrinos. Oscillations or decays of the heavier sterile neutrinos 
in the neutrino minimal standard model can give rise to the requisite lepton asymmetry, provided 
the masses of the heavier neutrinos are sufficiently degenerate. We study the renormalization group 
evolution of the mass difference of these singlet fermions to understand the degree of necessary 
fine-tuning. We construct an example of the model that can lead to a technically natural realization 
of this low-energy degeneracy. 
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I. INTRODUCTION 

The neutrino minimal standard model (z^MSM) is a 
renormalizable extension of the standard model (SM) by 
three singlet Majorana fermions with masses below the 
electroweak scale. As has been elaborated in Refs. [H-Il] 
(for a review see Q), this model enables us to solve, by 
economic means, four observational problems of the SM. 
Owing to the low-energy see-saw mechanism the i^MSM 
leads to nonzero masses of the active neutrino flavors and 
thus to neutrino oscillations. The model provides a can- 
didate for dark matter (DM) particle in form of a sterile 
neutrino in the mass range of 1 — 50 keV (see iS] for 
a review) . The coherent oscillations of the heavier neu- 
tral leptons insure the generation of baryon asymmetry 
of the Universe 0, [l^ • ^^'^ lastly, a (large) nonminimal 
coupling of the SM Higgs field with gravity leads to in- 
flation consistent with the cosmological observations [1] . 

In the ^MSM, the dark matter sterile neutrino is pro- 
duced at temperatures ^ 100 MeV due to the mixing 
with active leptonic flavors. The spectrum and the num- 
ber density of produced particles depends essentially on 
three parameters: the mixing angle with ordinary neu- 
trino, the sterile neutrino mass and the lepton asymmetry 
of the Universe (A£) at the production time (for a re- 
cent analysis, see [H], [13; earlier considerations can be 
found in [6-9]). The comparison of theoretical predic- 
tions with cosmological and astrophysical observations 
(such as Ly-a data and x-ray observations, constraining 
the free streaming length of DM particles and their mix- 
ing angle with neutrinos, respectively) lead to the conclu- 
sion that the (low-temperature) lepton asymmetry must 
be much larger than the baryon asymmetry AB ^ 10^^", 
^ > 3 X 105. 
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^ A possible way to evade this requirement is some modification 
of the I'MSM, aUowing the DM sterile neutrino interactions with 
other new particles, see, e.g. [l3- fl6| . Yet another possibility is 
related to primordial Higgs-inflation [T7| . 



As was shown in [3| , the presence of a pair of nearly de- 
generate heavier neutral leptons in the i^MSM may lead 
to production of the requisite large lepton asymmetry 
below the electroweak scale without a conflict with ob- 
served small baryon asymmetry (generated by the same 
particles and by sphalerons at electroweak temperatures) . 
Basically, the resonant production of AL occurs at de- 
coupling or during decays of singlet fermions which are 
taking place well below the sphaleron freeze-out. 



The requirement of generation of sufficient lepton 
asymmetry leads to the stringent conditions on the pa- 
rameters of the model, partially analyzed in Q. The 
most important of them is the level of the degeneracy 
of the pair of neutral leptons, which demands a severe 
fine-tuning of the masses and couplings of the z^MSM. 



The aim of this paper is to study the stability of nec- 
essary fine-tuning against radiative corrections (only the 
tree-level analysis has been made in Q). In particular, 
we will study the renormalization group (RG) evolution 
for the mass difference of the singlet fermions and formu- 
late the conditions that can lead to a technically natural 
realization of the low-energy degeneracy, required for the 
low-temperature resonant leptogenesis, essential for DM 
production. 



The paper is organized as follows: In Sec|TT]we review 
the basic structure of the i/MSM and its parametriza- 
tion. In Sec IIIII we explain the necessity of the degen- 
eracy between neutral leptons. In Sec IIVI we describe 
and analyze the renormalization group evolution of the 
essential parameters of the model. In Sec |V] we revisit 
different possible scenarios for singlet fermion mass split- 
ting and in Sec IVII we discuss an extention of the i^MSM 
by higher-dimensional operators. In Sec IVIII we discuss 
the phenomenological bounds obtained as a result of the 
fine-tuning. Finally, in Sec IVIIH we summarize our re- 
sults. 
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II. THE //MSM AND ITS PARAMETRIZATION 

We use the Lagrangian of the z^MSM in the following 
parametrization [3|, |l8| : 



^uMSM — Co + AC, 



(1) 



Co = CsM + J2 Niid^-r'^Ni - {K2La.N2^ 



1=2,3 



+ MN^Ns + h.c), 



(2) 



AC 



-hasLaNacj) - 



AM 



E 

7=2,3 



N}Ni + h.c, (3) 



where Nj are the right handed singlet leptons (/ — 2,3), 
(j) and La{a = e, /z, r) are the Higgs and the lepton dou- 
blets, respectively, /i is a matrix of the Yukawa coupling 
constants, M is the common mass of the two heavy neu- 
tral fermions, AM is the diagonal element of the Majo- 



rana mass matrix, and di — e 



M and AM are taken 



to be real. We have omitted the DM sterile neutrino Ni 
from the Lagrangian as its influence on the problem we 
are interested in is negligibly small (see Q for details). 

As has been shown in [3], one can solve for the active 
neutrino masses explicitly: 

2 

m=[0,^[F2F3±\h^h\23]}, (4) 

where = [h'^hlu and v ~ 174 GeV is the vac- 
uum expectation value of the Higgs field. One encoun- 
ters two different cases, namely the "normal hierarchy," 
mi = 0, 7712 = msoi ^ 0.009 eV, nis = matm — 0.05 eV, 
and the "inverted hierarchy," mi ~ m2 ~ matm — 
0.05 eV, ma = 0, |mi-m2| w ml^J {2mat7n) ^ 
8 X 10^^ eV. Normal hierarchy corresponds to the case 
when |ft.^/i|23 ~ F2F3, and the inverted hierarchy to the 
case when |ft.^ft.|23 ^ F2F3. Here mgoi, matm are the 
solar and atmospheric neutrino mass differences (for a 
review see jHI), 

(5) 
(6) 



= ''■65l°;23 X 10"^ eV^ 



-2 _ 2 40+0-12 



^atm 



IQ-^ eV^ 



For the lepton and baryon asymmetries of the Uni- 
verse, one of the essential parameters is the ratio 



^2 ^ (h^hhs 



(7) 



measuring the relative strength of the coupling of the 
neutral leptons A^2,3 to the active flavors. Without loss 
of generality it can always be chosen in the region e < f . 

The most important parameter is the mass difference 
between the two heaviest neutrinos. Successful baryoge- 
nesis and leptogenesis necessitates the mass difference to 



be small. At the tree level, there are two contributions 
to the mass difference: one related to the Majorana mass 
matrix and the other to the Higgs vacuum expectation 
value and Yukawa couplings Q : 



dmtre.e. = 



M 



where 



2{h'^h) 



23V 



2MAM 



(8) 



(9) 



Note that 2\h'^h\23v'^ = M\Amy\, where Am^ is the 
difference between active neutrino masses, and Am^, ~ 
0.04 eV (8 X fO"''^ eV) for normal (inverted) hierarchy 
(Am^ by definition is a physical quantity and is the RG 
invariant) . 

The analysis of radiative corrections to Eq. © and its 
RG evolution plays a central role in our considerations 
(see also accompanying paper [20j). 



III. THE FINE-TUNINGS NECESSARY FOR 
LOW-TEMPERATURE LEPTOGENESIS 

It has been shown in [3] that the singlet fermions iV2.3 
enter into thermal equilibrium at some temperature r+ ^ 
Mw {Mw being the intermediate vector boson mass) and 
freeze out later at T^>M. Then they decay at Td<M. 
The low-temperature lepton asymmetry can be generated 
at T ~ T_ and at T ~ T^j, to satisfy Sakharov out of 
equilibrium conditions. 

The Yukawa coupling constants in the j/MSM are very 
small, 



F2F3 — 



X 10 



-16, 



M 
GeV 



(10) 



where k — 1 (k — 2) for normal (inverted) hierarchy. 
Therefore, the production of substantial lepton asymme- 
try can only be possible in case of resonant N2 O N3 
transitions. This means that the oscillation rate Tqsc 
should be of the same order as the scattering rate 



^G^pT^matm 

7m 



(11) 



if the leptogenesis occurs at T ^ T_ , or the decay rate 



r 



N 



IQGlM^matm 

1927r3e 



(12) 



if T ^ Td. 

The quantity Fosc is related to the difference of physical 
masses of Majorana fermions dmphys as Tone ~ 5mphys, 
if T < M, or Tosc ~ MSmphys/T for T > M. So, to get 
a substantial lepton asymmetry at T ^ T_, we have to 
require that 



Mdmphys _ 
T- ~ Mo 



(13) 
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where Mq w Mpi/lM^/g^, Mpi = 1.2 x lO^^ GeV 
and the temperature dependence of the effective number 
of massless degrees of freedom g^/f may be taken from 
[iij. From ([13]), for T - M - 1 - 10 GeV we get, 

SmphyAT-) < 10^^^ - IQ-^^ GeV . (14) 

An even stronger condition must be true for asymmetry 
generation in N decays. Taking again AI ^ 1 — 10 GeV 
we have 

SrUphysiTd) < 10-23 - 10-19 GeV . (15) 

The physical mass difference of singlet fermions in Eqs. 
(|14I15P should be taken at the corresponding tempera- 
tures, r_ or Td- 

It is the smallness of these numbers in comparison 
with the observed neutrino mass difference Am,y (at least 
8 X 10~^^ GeV for the inverted hierarchy case) that leads 
to the fine-tuning problem. Indeed, there are several 
terms of different nature which make up the physical 
mass difference: the tree level Higgs contribution ^ Am,^, 
the tree level Majorana contribution ^ AM, the zero- 
temperature loop corrections, and finite-temperature cor- 
rections. The latter effects were analyzed in detail in 
3, Sec. 5.1. It was found there that these corrections 
lead to the induced finite-temperature mass difference 
Sniphys which is of the order of the Hubble rate H at 
temperature T_ (see Eq. (7.32) of ^) or smaller than 
H at Td- In other words, they do not spoil the resonant 
character of low-temperature leptogenesis, provided the 
zero-temperature mass difference is tuned to small values, 
given in Eqs. (|14I15I) above. Therefore, we will concen- 
trate on zero-temperature contributions in what follows. 
A part of the zero-temperature loop corrections can be 
absorbed into the tree-level Higgs contribution, convert- 
ing it to the physical mass difference of active neutrinos 
[Sj: 

Am M"^ 

5mphys = |Am,e^" -h AM| + 0{—^ — ) , (16) 

where a is some phase. Numerically, the last term in (fT6|) 
is of the order of 10"!*^ GeV (for M ~ 1 GeV), which 
is of the same order as ((Ti)) and is much larger than 
([15]). Therefore, the tuning of the tree- level Majorana 
contributions to the physical value Am^ must be done 
together with radiative corrections in order to achieve 
the required degeneracy Sniphys ^ Am^. 

Yet another attitude can be used in discussing the fine- 
tuning. Suppose that the compensation of the tree con- 
tributions in ([HI is associated with some symmetry which 
may potentially exist at some high energy scale, such as 
the Planck mass Mpi . Then the physical mass-difference 
SrUphys, which is determined by the low- energy parame- 
ters, is not in general zero, due to the running of these 
parameters. The consideration of the running of the rel- 
evant parameters, together with computation of the ra- 
diative corrections [20| , would allow then to estimate the 
"natural" values of the mass difference. This is the pur- 
pose of the next section. 



IV. RG EVOLUTION OF THE 
MASS-DIFFERENCE 

The RG running of the Majorana masses and of sterile- 
active Yukawa couplings can be extracted from [2l|, 

{Aixfj^MR = {h^h)MR + Mnih^hf , (17) 

{Airf^h = {hh^ - IyJY, - ^gl - -al 
^ ' dt 1 2 2 4 1 4 2 

+ Tr[3yjr„ + WlYd + hh^ + Y^Y^i^h, (18) 

where i = In ^ (for definiteness we take /xq = top quark 
mass), Mp is the Majorana mass matrix, Yg is the diago- 
nal charged lepton matrix, Yu(d) is the Yukawa coupling 
matrix for the up (down) quarks, and 51, g2 represent the 
gauge coupling for U(l) and SU(2), respectively. Since 
^ contains the vacuum expectation value of the Higgs 
field, v'^ = rnjj/2X {m?jj is the mass parameter in the SM, 
A is the scalar self-coupling)0, we will also need the RG 
runnings of m?^ and A. They are given, for example, in 




6y?), (19) 



(47r)2^ = (12A-35?-95|)A-|-^3^ 

+ \9l9l + \9t + l'2{\-yl)vl (20) 

To complete the system, we add the RG equation for 
the top quark Yukawa coupling yt , 

(""'f = »'(5=''-i''-i'-«»»)- <^'> 

The RG evolution of other parameters of the SM will not 
be needed in this section. 

Using the fact that h <C 51,2, and neglecting numer- 
ically small Yukawa couplings of charged leptons and 
quarks (except j/j), we get the RG equation for the Ma- 
jorana mass difference and for h: 

(47r)2^AA^ = 2M|/it/i|23 , (22) 
{A.f±h = h{?.yl-\gl-\9l]. (23) 



^ In this paper, we use the Higgs potential form as: 

y(</.) = -/.2(<At</.) + A(^t^)2. 
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We consider the Yukawa couplings for the sterile neutri- 
nos to be real in this section. This simplifies the com- 
putations, but does not change the numerical estimates 
and qualitative conclusions. 

So, the change in AAf due to RG evolution from the 
Planck scale to fio is of the order of 

^"^''Inr^^ ^.(6.5x10-- -4x10--)-^' 



(47r)2 



Mo 



GeV ' 
(24) 

depending on the type of neutrino hierarchy. This num- 
ber for the inverted hierarchy is roughly of the same order 
as (fT4| . but exceeds considerably (fT5|) . 

The RG evolution of the Higgs contribution to the mass 
difference is much more substantial. In order to make an 
estimate of this contribution, we fix the initial condition 
for the Yukawa couplings at the Planck scale and run 
them down to the electroweak scale. For definiteness we 
take the inverted hierarchy case, the conclusions for the 
normal hierarchy are the same. 

It was shown in |^ that for inverted hierarchy the cou- 
pling of A^2 to T flavour (consideration of other leptonic 
families do not change the results) is 



cosyi2e 



I sin6'i2e 



(25) 



where tan^ 6112 — 0.48 from active neutrino oscillation 
data [1^ and C is a phase factor that varies between 
and 27r. A little algebra shows that 



0.14, 



hr2 



F2 



0.7 



(26) 



Similar relations hold for |/it-3| with F2 replaced by F3. 

Using and ([TU]), it follows that the maximum and 
minimum values of the Yukawa couplings, for AI — 1 
GeV and e = 1, are 



3.4 X 10" 



\hT2 



1.7 X 10" 



(27) 



A similar set of maximum and minimum values hold for 
\hr3\. The values of Yukawa couplings scale as ha2 oc 



y/M/e, ha3 oc 

Using the above-mentioned bounds along with Eqs. 
(HH),® and ([23]), we obtain the RG evolution of the 
Higgs contribution down to the electroweak scale. Per- 
forming the necessary numerical computation, we find 
that the Higgs contribution to the mass-difference ranges 
from 10^1° GeV to 10"" GeV depending on whether the 
maximum or minimum value is chosen for the Yukawa 
couplings. This value is of the order of Am^, meaning 
that if in (jH]) is tuned to zero (implying Smtree = 0) 
at the Planck scale, the physical mass difference between 
N2 and N3 will generically be of the order of Am,y, ex- 
ceeding considerably the required values (I14ll5p . A way 
out is a fine-tuning of the Higgs mass, which can be cho- 
sen in such a way that Smtree{Mpi) = Smtree{Mw) = 0. 
In Fig. [1] we give the RG evolution of the mass-difference 
drritreeil^) for a particular choice of Higgs mass (145 
GeV), which leads to this situation. 
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FIG. 1. RG evolution of Smtree for a Higgs mass of 
145 GeV. 



V. SCENARIOS FOR SINGLET FERMION 
MASS DIFFERENCE REVISITED 

To summarize, if there were a compensation in the 
tree-level contribution to the mass-difference dnitree due 
to existence of some symmetry at a higher-energy scale 
(Mpi), it does not exist at the electroweak scale and the 
physical mass difference acquires a nonzero value, which 
is of the same order as that of the change in Higgs contri- 
bution: 10~i° - 10"" GeV. The reason is that the Ma- 
jorana contribution to runs much more slowly than 
the Higgs one. One can arrive to (|14I15I) either by tuning 
of the Higgs mass, or by tuning the initial condition for 
AM at the Planck scale. In Sec. IVII we will discuss how 
this situation can be changed in some extension of the 
j/MSM. Meanwhile, we will describe the possible scenar- 
ios (sf. [3]) for singlet fermion mass difference account- 
ing for RG behavior studied in Sec. |lVl Depending on 
the relative importance of the different contributions to 
SiTiphys, they may be classified as follows: 

• Scenario la: AM{Mpi) — 0. In this case the phys- 
ical mass difference is mostly due to the tree-level 
Higgs condensate and loop corrections. One can 
easily check from Eqs. (|H), ^ that this leads to 



SiTlpf^yg 



Am,, 



Thus, SrUphys — matm - m^oi 
normal hierarchy and Smriphyg 
GeV for inverted. 



(28) 

5 X 10~" GeV for 
8x10-13 



• Scenario lb: This corresponds to the situation in 
which the tree-level Higgs and the Majorana con- 
tribution are of the same order of magnitude at 
the Planck scale, including the case when there is a 
compensation of the two contributions, rin?{Mpi) — 
0. In this case and without any special fine-tuning 
one gets 



Am,. 



(29) 



Scenario II: The physical mass difference of the 
singlet fermions is much smaller than the active 
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neutrino mass difference, i.e. 



(30) 



Only this scenario can lead to production of sub- 
stantial low-temperature Icpton asymmetry and 
thus to resonant sterile neutrino dark matter pro- 
duction. It requires a fine-tuning between contri- 
butions of the different nature (Higgs, Majorana, 
and loop corrections) and thus is "unnatural" in a 
technical sense 0. 

• Scenario III: The Majorana contribution domi- 
nates, in which case 



Srriphys > Amy. 



(31) 



The scenarios la, lb, and III are natural in the techni- 
cal sense. However, they do not lead to low-temperature 
lepton asymmetry. In Sec. IVII we will present an exten- 
sion of the i^MSM in which Scenario II can be realized 
as a natural possibility. 



VI. HIGHER-DIMENSIONAL OPERATORS 

The z/MSM can be extended in several ways, which 
may break the relations (|4|8|9p . being the basis of the 
analysis of the previous sections. The simplest possibility, 
in the spirit of effective field theories, is to add higher- 
dimensional operators. There are five independent five- 
dimensional operators which can be constructed from the 
fields of the z^MSM. One contains the fields of the SM 
only. 



(32) 



The other two include the singlet fermions as follows: 



ci'N'iNj4,'^(^, c'/id^Nj^d^Nj. (33) 
Yet another two include L and N simultaneously. 



These initial conditions require that the active neutrino 
mass hierarchy must be inverted, as follows from Eq. Q . 
Then, if charged lepton Yukawa couplings Yg are set to 
zero, the Lagrangian of the z/MSM has an extra global 
leptonic U(l) symmetry (ISj]. This means, that the re- 
lations ([35]) are RG scale-independent, and that the sin- 
glet fermions will remain exactly degenerate. In fact, 
the charged leptonic Yukawa couplings violate explicitly 
this symmetry and thus lead to the breaking of exact 
degeneracy at small energies. So, some mass difference 
between the singlet fermions will be generated. As we 
will see below, it is much smaller than the observed mass 
difference between active neutrinos. The inclusion of the 
operator (f32|) is needed to generate the observed active 
neutrino mass difference and requires ci ~ l/(10^^Gey); 
the common mass of active neutrinos is due to -/V2.3- 

Let us make an estimate of the splitting between A^2,3 
owing to charged lepton Yukawas. Among Ye the tau- 
lepton coupling is the largest, and Y^Ye, which was pre- 
viously omitted from Eq. p^ . has the form 



Y^Y, 



^0 

,0 y2 



(36) 



The evolution of AAI and (/i^/i)23, following from Eqs. 
(|17I18I) with initial conditions (I55|) is shown in Figs. 
and [3] respectively. We give only the plots correspond- 
ing to the minimum value of the singlet fermion Yukawa 
couplings [see ^7} ] and for M = 1 GeV. The energy 



3. X 10- 
2. X 10- 
1.x 10- 




(34) 



Here 
GeV 



Ci are new coupling constants with dimension 
^. The operator (|32p can change the relation (U), 
whereas operators (|33l34p can change (|8l9p . 

There are several possible ways how the Scenario II 
can be made technically natural. The first one is based 
on the use of the operator ([5^ . which contributes to the 
active fermion mass difference. The idea is as follows. 
Suppose that at Mpi the following two relations hold si- 
multaneously: 



AM = , (/i^/i)23 - 



(35) 



^ By "technically" natural we mean the situation in which the fine- 
tuning made at the high energy scale persists to small energies. 



FIG. 2. RG evolution of AM. 

scale has been varied from the top quark mass (t — 0) 
to Planck scale {t = 39.14). We can see that both \AM\ 
and {h^h)23 decrease steadily with the increase of t. 

Numerically, for the minimal value of /i-r2 [case (a) for 
future reference], the mass difference is given by (sf eq. 



Sm 



phys 



5x 10-^V"-1.2x 10"" 



GeV , (37) 



GeV^ 

establishing a typical scale ^ 10~^^ GeV. 

Analogous computation for the maximal value of hr2 
[case (b)] leads us to the similar dependence: the mass- 
difference is given by 

Af2 



1.3x 10-^V"-3x 10" 



16 



GeV 



GeV, (38) 
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FIG. 3. RG evolution of \li^h\23. 



with a typical scale of ^ 10^^® GeV. 

In order to find the different possible values of mass- 
difference, we analyze the two cases separately for differ- 
ent values of M. 

• case (a): From Eq. ([57)) we see that for a partic- 
ular choice of the Majorana mass (M ~ 2 GeV), 
the mass difference can be made arbitrarily small 
by a suitable choice of the Majorana phase (a). In 
particular, the mass difference vanishes for a = 0. 
Therefore, it may be safely concluded that for 
M ~ 2 GeV, the following bounds hold true for 



< 5mphys < 10"^^ GcV. 



(39) 



We now turn to the two other possibilities: M ^ 2 
GeV and M > 2 GcV. Considering M < 2 GeV 
first, we can see that 



Smphys = 5 X 10 GeV, 
while for M > 2 GeV, 



6mphys = 1-2 X 10 



-17 



GeV 



(40) 



(41) 



• Case (b): Performing similar analysis, we see that 
once again for M ~ 2 GeV, the mass difference be- 
comes arbitrarily small, depending on the choice of 
the Majorana phase. We thus establish the bound 
for Smphys in this case: 



< Smphys < 2.6 X 10"^^ GeV 
For M< 2 GeV, 

Smphys = 1.3 X IQ-is GeV, 
while for M > 2 GeV, 



Sm^ 



■phys 



= 3 X 10 



-16 



GeV 



(42) 



(43) 



(44) 



Thus, we see that the mass difference Smphys so ob- 
tained in the two cases is indeed much smaller than the 
active neutrino mass difference, what is needed for the 



low-temperature resonant leptogenesis. As we have al- 
ready said, the role of the higher-dimensional operator 
([32l) is to provide the additional contribution to mass dif- 
ference of the active flavors. The amplitude of all other 
operators must be small enough in order not to spoil the 
relations p7|38|) . which is technically natural. In Sec. 
IVIII we will consider different bounds on parameters of 
the i^MSM which appear in this scenario. 

Yet another possible way to make the Scenario II 
technically natural is based on operators (|34l) . These op- 
erators contribute to the singlet fermion mass with the 
terms of the order c^^hajv'^, c^^hajv"^. If the RG run- 
ning of some combination of these contributions happens 
to be the same as the running of the Higgs contribution 
2{h^h)23v'^ /M, this can be used to cancel the largest ef- 
fect ([24]) . If true, the physical mass difference can be 
"naturally" made of the order of the last term in Eq. 
P^ . The analysis of this possibility goes beyond the 
scope of the present paper. 



VII. PHENOMENOLOGICAL BOUNDS 
PARAMETERS OF i^MSM FROM 
NATURALNESS 



ON 



Let us assume that the physical mass difference is in- 
deed given by the relations derived in Sec. |VT1 Is the 
baryogenesis due to singlet fermion oscillations still op- 
erational? What can be said about the parameters e and 
M from the requirement of resonant dark matter produc- 
tion? 

To answer the first question, let us determine the cru- 
cial parameter of the baryogenesis - the number of os- 
cillations x(Tsph) of singlet fermions before the freezing 
of sphalerons at temperature Tgph ~ 150 GeV The 
high-temperature mass difference of the singlet fermions 
comes from two-loop graphs which include the square of 
charged lepton Yukawa coupling and square of hai and 
corresponds to the double Higgs exchange. It can be es- 
timated as 



leading to 



AM{T)^ 



hr2VT 

16 



16 



sph 



(45) 



(46) 



Yet another contribution to AM(T)^ comes from the 
Higgs condensate. It is of the same order of magnitude as 
(j45l) . since at the sphaleron freeze out v{T) ~ T. Numer- 
ically, a;(Tsph) ~ (10"'' — 10~^)/e. Since baryon asym- 
metry is proportional to x for x 1, it is automatically 
smaller than the low-temperature lepton asymmetry by 
a factor ^ 10^ (for e ^ 1), providing a potential explana- 
tion of this hierarchy. At the same time, with this value of 
a;, the baryon asymmetry can be as large as AB ^ 10~^, 
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O.ISk-B 

a;(T_) = — {GFMofniatmSmphysiT) 



exceeding the observed one. In other words, the answer x{T^)/2'k of iV2,3 given by 
to the first question is positive. 

To analyze the second question, we wih consider two 
choices of ft.1-2, referred to as (case (a) and case (b)) 
above. As has been shown in the lepton asymmetry 
can be generated in either by decays or oscillations of the 
sterile neutrinos. We will consider these two situations 
separately. 



(51) 



A. Leptogenesis by decay of sterile neutrinos 

Considering the case of decays first, it has been shown 
in Q that the maximal lepton asymmetry A which can 
be generated in decays of A^2,3 is of the order 



MoSmphys 



(47) 



where A^ax = 2/11. The condition that the decays of 
A^2.3 occur above the dark matter creation temperature 
~ ioO MeV reads 



M 
GeV 



> 1.4 



2 X 10-3 



(48) 



To produce a necessary amount of dark matter sterile 
neutrinos, it is required that (l^ 



A > 2 X 10' 



(49) 



• Case (a): When M < 2 GeV, using (gOl) , (gll) and 
(j48|) . we arrive at the minimum value of M to be 
3 GeV, which contradicts the starting assumption 
that M < 2 GeV. On the other hand, for M » 
2 GeV, using (j4T1) and (|47l) . we arrive at a lower 
bound on e: 



(50) 
is 



e>0.1. 

Corresponding minimum value of M [from 
- 4 GeV. 



• Case (b): Performing a similar analysis as before, 
we once again reject the case when M <^ 2 GeV. 
For the case when M ^ 2 GeV, we obtain a lower 
bound on e to be 2.3, conflicting with the condition 
e < 1. 

To summarize, a sufficient lepton asymmetry in decays 
of iV2,3 can be generated in case (a) for M>4 GeV and 
e>0.1. Or, one has to require the fine-tuning AI ~ 2 
GeV and a ^ 1, making the physical mass difference 
even smaller than ~ 10~^^ GeV. 



B. Leptogenesis from coherent oscillations of 
sterile neutrinos 

As has been shown in Q , the most important parame- 
ter which determines the value of lepton asymmetry gen- 
erated at temperature r_ is the number of oscillations 



where k is 1 or 2 depending on normal or inverted hi- 
erarchy, _B = 5 and Gf is the Fermi coupling constant. 
The case SniphysiO) = 0, in which 5mphys{T-) 7^ due to 
finite-temperature effects, leads to x{T-) = a;^ — 10 and 
to some region in the e, M plane leading to the required 
lepton asymmetry (j49p . This analysis stays in force also 
if 6mphys{0) 7^ 0, provided 



Smphys{0) < 



0.l5KB{GFMQ)^matm 



(52) 



We will give below the corresponding parts of the phase 
space. For this end we find x(T_) from Eq. (j51l) replacing 
Smphys{T) by drriphys defined in Eqs. (|39II44I) . 

• Case (a): For M < 2 GeV, using PO ]) and 
we get 

0.25 



x(T_) 



(53) 



which implies e > 0.025 to satisfy inequality 
x{T^) < XT- For M > 2 GeV, using ^ and 
dSlI), we get 



x(r_) 



0.06 



(54) 



e GeV^ 

which implies e > O.OOGAf^/GeV^. Plotting e as a 
function of M, for the two situations, the allowed 
values of e — M are depicted by the shaded portion 
in Fig. m 




M 

10 GeV 



FIG. 4. Shaded region of the plot represents the values of e-M 
in which the analysis of Ref. 0| is not changed for the case 
(a). The e axis is shown in logarithmic scale. 

• Case (b): Similar analysis for M ^ 2 GeV leads 
to 



x{T-) = — , 
e 



(55) 



which implies e > 0.64. Again for m ^ 2 GeV, we 
arrive at the relation: 



x(T_) 



1.48 M2 
e GeV^ 



(56) 
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which imphes e > 0.148Af^/GeV^. These regions 
are shown in Fig. [S] 



1. 

0.95 
0.9 

0.85 
0.8 

0.75 

0.7 

0.65 
0. 



M 
4.0 GeV 



FIG. 5. The same as in Fig. 4 for the case (b). In this plot, 
we only show the region where e < 1, which we assume to be 
true throughout the paper. 



As in Sec. IVII A| the case when M ~ 2 GeV and a < 1 
is special, and all values for e and Sniphys, found in Sec. 
IVII are allowed. 



tion. In this paper we studied the RG evolution of the 
mass splitting between neutral leptons iV2,3 of the i/MSM, 
essential for resonant production of low-temperature lep- 
ton asymmetry, followed by resonant production of dark 
matter sterile neutrinos. We found that the mass dif- 
ferences of the order or greater than the observed mass 
differences in the active neutrino sector are natural in 
the technical sense. In other words, the RG running of it 
from the Planck to the low-energy scale leads to correc- 
tions ~ Amy. 

At the same time, the low-temperature resonant lep- 
togenesis requires the splitting that is much smaller than 
^ Ami, and thus is "fine-tuned." We described an ex- 
tension of the z^MSM by higher-dimensional operators, 
in which this fine-tuning is due to an approximate sym- 
metry of the theory at the Planck scale and which is not 
spoiled by the RG evolution. It requires the hierarchy 
of neutrino masses to be inverted. We analyzed the con- 
straints on the masses and couplings of the singlet leptons 
in this scenario and demonstrated its feasibility. 



VIII. CONCLUSIONS 
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